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ABSTRACT 

Theory for aeolian tone intensity created by rigid 
circular cylinders has been found to agree quite veil with 
experimental results. This investigation is to experimentally 
explore the wake of a cylinder shedding vertices at a frequency 
that is equal to one of the natural vibrational frequencies of 
the cylinder, Inves tigation v;as carried out in the following 



areas : 






t 


i. 


lateral spacial correlation of vortex 
shedding. 




2. 


distance behind cylinder and value of 
RMS velocity at assumed position where 
vortices start to n roll up"* 




3. 


observation of the phase between velocity 
signals at two points along a cylinder. 
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CHAPTER I 



INTRODUCTION 



The study of fluid flow past circular cylinders has been 
in progress since the latter part of the 19th century when Strouhal 
first showed that the vibrations causing the sound were transverse 
to the direction of the wind. The early interest in the phenomena 
was mainly concerned with the vortex shedding frequency as a 
function of Reynolds number. 

Until recently, very little work has been done concerning 

the intensity of vortex shedding sounds known as aeolian tones, 

Phillips ^ and Etkin, Korbacher and Keefe ^ have derived 

theoretical expressions for the radiated sound intensity by starting 
(3) 

with Lighthills basic theory of aerodynamic sound as had already 

(/,) 

been extended to rigid bodies by Curie , 

Phillips used the product of Z Z^ to account for the 
fact that the force may not be coherent over the entire length of 
the cylinder. The use of the product Z Z^ is valid if the body is 
long compared to the coherent length. His resulting equation for 
sound intensity is: 

1 « S 2 C 2 It 3 U 3 l t sin 2 0 

16 r 2 * 

( 1 ) 



1 



(-) “ indicated time averages 

p « free stream density 

r,0 ~ polar coordinates (9 being the upstream direction), 

S w = Strouhal number 

1 o 

d Q cylinder diameter 

N “ vortex shedding frequency 

U - free stream velocity 

K lift coefficient 
M ~ Mach number 

l - length of cylinder shedding vortices 

- spanwise coherent length of vortices 
I - acoustic intensity 

Etkin, Korbacher, and Keefe derived a more general equation that 

considers any length of coherent force. 

In order to allow for longer lengths of coherent force , 

the centroid of the area under the correlation coefficient must 

2 

be included. The effective length that is substituted for £ is 

l (&~y) rather than £ £ , 
c c 

Y w centroid of the area under the cross correlation 
coefficient curve, y becomes small as the length of coherent 
force becomes small. Therefore, equation (1) is apporached for 



2 



small A , The more general equation from Etkin, et al* shows: 



I 



.JL,- 

16 2 

r 



M 3 U 3 £ (£~y) 



sin^B 



( 2 ) 



Experimental investigation of radiated sound intensity 
from circular cylinders when shedding vortex wakes has been made 
by Gerrard^ 3 ^ and Keefe In general* these experimental re- 
sults confirm the theory for rigid bodies especially in the 
lower Reynolds number region* 

The flow past cylinders that are non-rigid bodies has 
not been investigated sufficiently to determine what happens when 
the vortex shedding frequency and a natural frequency of the 
vortex producing body coincide* This thesis is a report of 
experimental investigation of certain properties of the wake 
behind a cylinder shedding vortices at a natural vibrational 
frequency of the cylinder producing the vortices* 

The theoretical equation mentioned earlier for sound 
intensity can be simplified for evaluation to: 



1 « c 2 (*-YH 

* C 



( 3 ) 



which shows the three terms that may change for the situation 

of frequency coincidence C g , and Y* Current work using the 

~ c 



3 



same experimental set up is proceeding to obtain experimental values 
for sound intensity, lift coefficient > and any "locking together" 
of the two frequencies r Phillips ^ expresses the opinion that the 
correlation length may increase greatly when the vortex shedding 
frequency and the vibration frequency coincide due to a coupling 
of fluid motions at different axial positions caused by the cylinder 
motion. This supposition was explored experimentally and the re- 



sults are presented here 



CHAPTER II 



PROCEDURE 



1 1 . 1 , r £^^A° l \ 



The aeolian tone test rig, which is pictured in 



Appendix A, was used to hold a cylinder of *154 inches diameter 
and 15 inches in length under very high tension to produce the 
desired natural frequency* It was desirable to have one of the 
lower harmonic vibrational frequency coincide with the vortex 



cylinder's vibration as possible* The cylinder was put under as 
high tension as the test rig was physically capable of holding. 
The table below gives the natural frequencies of vibration and 
the air speed and Reynolds number required to cause vortices to 
shed at these frequencies. Only the wake of the 2nd, and 3rd, 
harmonic was explored. 



shedding frequencies in order to eliminate as many nodes in the 



N (cps) U (fps) 




2nd, harmonic 



62,5 



4,57 x 10 



3 



1030 



3rd, harmonic 1520 



83 



6,07 x 10 



3 



5 



The cross-correlation coefficient was computed in the 
following manner: 




( 4 ) 

R (z) 4:3 the cross-correlation coefficient as a function of 
spacial separation of the hot wire anemometers. 

The cross-cor relation coefficient is stated here to be 
only a function of the separation and not a function of where the 
separation is located along the cylinder. The investigation to 
determine if this coefficient and, in turn, the correlation length 
is independent of location along the cylinder was not carried out 
and is a subject for future investigation* 

The values of R(z) were then plotted against z f See 
Figures [4], [5], [6], and [7]* A mechanical integrator was used 
to determine the area under the R ( z ) curve. This value is defined 
as one half the correlation length* 

i.c. * c » j " K a:j (,) dz 

•~co 

(5) 

A plot of a? a function of Reynolds number is shown on Figure [8], 
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In this range of Reynolds numbers , the Strouhal number is practl- 
cally constant at *205, The cylinder consisted of a ,091” diameter 
steel wire over which was placed a Teflon tube that brought the 
overall diameter to .154 inches t 

This cylinder was placed in the open jet test section 
of the low turbulence wind tunnel in the M,I,T, Acoustics and 
Vibrations Laboratory, Figure [1A] shows the dimensions of the 
wind tunnel. It is more fully described in a report by Hanson 
For the correlation measurements one of the hot wire 
sensing elements was placed in a stationary position while the 
other sensing element was held on a traverse mechanism that could 
be moved in a direction parallel to the cylinder. The two hot 
wires were oriented parallel to the cylinder and positioned at 
the edge of the wake where the signal is periodic f but not 
turbulent. This position varied from speed to speed F but in 
general was about 1/2 inch back and 1/4 inch to one side of the 
center of the cylinder* See Figure [2], The signals from these 
two Constant Temperature Anemometers (CTA) went to a correlation 
unit that produced a sum or difference output. The RMS output was 
recorded on a Graphic Level Recorder, Four outputs were recorded 

, V? V? 

See Figure (3 ]* A, !> = voltage signals from the two C'J'A's, 




To further check the correlation at coincidence of 



frequencies a dual beam oscilloscope was used to observe the 
phase relation of the signals from each hot wire* 

To ensure that the vortex shedding frequency is 
coincident with a natural frequency* a sine wave signal from an 
oscillator* set on the desired natural frequency by an electronic 
counter* was displayed on the dual beam scope* The speed of the 
tunnel was then adjusted to cause the hot wire signal to match 
that frequency of the oscillator* 

After setting the wind tunnel speed* the signals from 
both CTA’s were displayed on the oscilloscope for observation of 
the phase between the two signals. 
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1 1 , 2 » m , J fo x:uy ( m IMIS Vejx'o i t^ 

To determine the position where it is assumed that the 
vortices start to 'toll up" one hot wire anemometer was used with a 
two directional traverse mechanism. A search technique was used in 
the following manner. Starting very close behind the cylinder * the 
wake was traversed noting the position at which the maximum RMS 
velocity occurred. The value of RMS voltage and the DC voltage 
output of the anemometer was recorded at this location along with 
the position. The hot wire sensing device was then backed away 
from the wire a little further, The same traverse was made 
recording the same values at this point. This was repeated until 
the maximum value of RMS voltage was bracketed, Then> the exact 
position was located by searching within the bracketed space, A 
typical plot of how the maximum RMS voltage varies as a function 
of distance behind the wire is presented as Figure [9], The maxi- 
mum value is presumed to be the location of where the vortices 
start "rolling up", 

A traverse behind the cylinder was performed at three 
speeds to have an idea of the properties of the wake. The results 
of RMS velocity and mean velocity are plotted and presented as 
Figure [10 ], [11], and [12], 
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CHAPTER III 



DISCUSSION OF RESULTS 



1 1 1 . 1 , Lateral Spacl^ 

As stated in the introduction, the equation for the intensity 
of aeolian tones is proportional to lift coefficient, correlation length* 
and centroid of the area under the cross-correlation coefficient curve 
as thus? 



■Whwe* 

I - C £ 2 (Z-y)Z c 



( 3 ) 



The correlation length is the effective length of a 
coherent vortex* As the Reynolds number increases, the two 
dimensional qualities of the vortex shedding start to disappear* 

In other words, there is a span-wise variation in the phase of the 
vortex shedding* The average coherent length of shed vortices is 
used in calculating the aeolian tone intensity. This average length 
is obtained by a procedure using two-point correlation of velocity 
along a cylinder as described in Procedure* 

Correlation of velocity variations rather than force were 
made for two reasons* First, it was impractical to put pressure 
taps on the cylinder due to the small diameter* Secondly, because 
the cylinder has an induced motion the pressure tap readings would 
consist of a component due to the motion as well as the pressure 
caused by flow fluctuations* 



10 



el Baroudi 



made similar cross-correlation measurements on 



( 8 ) 

a rigid cylinder* His data ranged in Reynolds number from about 
4 4 

10 to 4,5 x 10 , He found that the correlation length increased 
as Reynolds number increased. His value of increasing from 
3, 2d at Reynolds number 104 to 6d at Reynolds number 4,5 x 10^ , 

Figure [8] t shows correlation length decreasing with 

increasing Reynolds number from about 8d at Reynolds number of 

3 3 

6,5 x 10 to 5d at Reynolds number of 11 x 10 , 

No apparent change in the trend of was noted at points 

of coincidence of the vortex shedding frequency and a natural vibra- 
tional frequency of the cylinder. It was thought that with the 

cylinder vibrating primarily in one mode that the coherent length 

l ( 9 ) 

would increase possibly approaching ^ (n 53 mode number), Gerrard , 

in fact* measured the correlation coefficient at low Reynolds number 

and found it to be of the following shape: 
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In the measurements obtained during this inves tigation, 
the value of the correlation coefficient approached zero monotoni- 
cally at all velocities including those at which the vortex shedding 
frequencies coincided with a natural frequency of the cylinder* 

More cross-correlation measurements should be made with 
particular emphasis on having the vortex shedding frequency exactly 
match the natural frequency and also more measurements at frequen- 
cies near the natural frequency* 

When observing the signals from both hot wires simultan- 
eously on a dual beam oscilloscope f with the hot wire probes 
adjacent to one another the signals remained inphase as would be 
expected* As the probes were moved apart parallel to the cylinder, 
the two signals did not remain inphase nor show a constant phase 
shift* Instead the two signals would move in and cut of phase with 
no apparent pattern* 

To check for phase reversal when crossing nodal points 
the hot wires were placed on each side of a nodal point at a point 
midway between the nodal points as thus for the 3rd, harmonic* 



CTA 
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Again the signal moved continuously in and out of phase 
as previously described * 

In all these observations^ on the dual beam scope* there 
didn ! t appear to be any change of observed signal whether the 
vortex shedding frequency was in coincidence or not with a 
natural frequency * 
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Ill .2. Position of Maximum RMS Velocity 



The position of the maximum RMS velocity is assumed to be 
the point where the vortex is "rolling up". 

The maximum RMS voltage referred to is the unfiltered hot 
wire anemometer signal. The signal at the position of maximum RMS 
velocity contained random frequencies together with the vortex 
shedding frequency, A frequency spectrum was taken of the hot 
wire anemometer signal. The background noise was at least 15db 
lower at all times than the signal at the vortex shedding frequency, 
Therefore* the addition to the broad band signal was less than 
0, 5db or ,01 volts,. 

As shown in Figure [9] the magnitude of the fluctuating 
velocity builds up linearly as the hot wire is moved downstream 
to a maximum then falls off showing the viscous diffusion of the 
vor ticity « 

(9) 

Birkhoff assumes that the wake acts as an oscillating 
flat plate attached to the cylinder providing alternating lift as 
the wake swings from side to side. He develops a satisfactory 
Strouhal number using this assumption. If it is assumed that the 
point of maximum RMS voltage is the point where the vortex starts 
"rolling up" and that Birkhoff f s evaluation of the wake acting as 
an oscillator providing alternating lift, then the distance 
downstream to the point of maximum RMS voltage and also the 



. 

I 



distance off to one side of centerline where the maximum PvMS voltage 
occurs would affect the lift produced by the wake* 

The distance downstream to the point of maximum RMS 
voltage is seen from Figure [13] to decrease linearly in the 
Reynolds number range considered* This follows the same trend as 
reported in Goldstein^"^ from data taken by Schiller and Linke* 

Again we see no deviation from this trend when the vortex 
shedding frequency and a nautral frequency coincide. 

Figures [10* [11 ]* and [12] show the mean velocity and 
RMS velocity of the wake at three different speeds with one of the 
situations having the vortex shedding frequency coincide with the 
3rd, harmonic of the natural frequency. The results show that 
with Reynolds number increasing the maximum difference in mean 
velocity increases indicating an increase in vortex strength. 

The maximum RMS velocity occurs at the point where the 
mean velocity is halfway between the maximum and minimum mean 
velocity. This is expected because this is the area of maximum 
shear flow. The slope of the mean velocity increases in the 
region of maximum shear flow with an increase in Reynolds number* 

In conclusion* there did not appear to be any abnormal- 
ties in the properties of the wake when coincidence of frequencies 
occurred * 
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APPENDIX A 



FIGURES 
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Test Cylinder 




FIGURE 1 

AEOLIAN TOME TEST RIG 
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TEST CYLINDER 



FIGURE 2 

SKETCH OF TliE LOCATION OF THE HOT WIRE PROLES IN 
RELATION TO THE TEST CYLINDER 



PROBE A 



PROBE B 
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CYLINDER 




EXPERIMENTAL SETUP FOR CORRELATION CURVE DATA 
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SAMPLE OUTPUT FROM GRAPHIC LEV El, RECORDER 
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ateral Spacial Correlation of 
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FIGURE 4 



Lateral Spacial Correlation of 
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Lateral Spacial Correlate 
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FIGURE 6 



Components of Cylinder Wake 
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MAXIMUM RMS VELOCITY AS A FUNCTION OF DISTANCE DOWNSTREAM 
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FIGURE 10 

RMS VELOCITY AND .MEAN VELOCITY AS A FUNCTION OF y 
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RMS VELOCITY ASD IliS/Vv VELOCITY AS A FUACTIOM Or y 
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FIGURE .1.2 

RMS VELOCITY AND IEAN VELOCITY AS A FUNCTION OF y 



29 



downstream from cylinder 
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DISTANCE TO "ROLLING U?" 0? VORTEX AND KAXII'TJM RMS VOLTAGE 
AS A FUNCTION OF REYNOLDS NUMLER 



APPENDIX T> 



EQUIPMENT USED FOR CORRELATION LENGTH MEASUREMENTS 
Thermo-Systems % Inc, 

Model 1010A Cons tant Temperature Anemometer 

Model 1015B Correlator 

Model 1031-2 Monitor and Power Supply with 
(2) Model 1035 Linearizer nodules. 

The sensing device used was a Tungsten Hot Wire , 00015 " 
diameter and ,05" long. 

The output of the Constant Temperature Anemometer was 
led into the linearizer modules so that given fluctuations in 
velocity would give the same voltage fluctuations for any mean 
velocity. The linearized output was sent to the correlator unit 
that summed and differenced the two signals; also the separate 
signals were obtained from the output of this unit. The output 
from the correlator was sent to the BxUel and Kjacr Frequency 
Analyser connected to the Graphic Level Recorder on vrhich the 
signals were recorded in db level. See Figure [3], The output 
was continuously observed on an oscilloscope. This was converted 
to RMS voltage and used in the cross-correlation coefficient 
equation. The magnitude of fluctuating velocity was not computed 
as the conversion factor cancels out of the R(z) equation. 
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The free stream velocity of the wind tunnel was obtained 
by placing the hot wire sensing device in a free stream portion 
of the open test section and recording the output voltage from the 
linear izer to use on the previously obtained calibration curve. 
From this velocity* together with the cylinder diameter and 
kinematic velocity* the Reynolds number was computed* 

The method used to check for coincidence of frequencies 
was to run a frequency analysis of the output of an accelerometer 
mounted on a support of the wire. When the two frequencies 
coincided* the natural frequency would be reenforced to its 
highest value* A more accurate method of matching frequencies 
was devised described on page (8) , 



32 



APPENDIX C 



RECOMMENDATIONS 

l t Take correlation coefficient measurements at different locations 
along the cylinder to determine if there is spacial homogeniety, 
2* Work with a system that will allow a greater amplitude of 
cylinder vibration to determine how amplitude of vibration 
affects the correlation length* The correlation lengths should 
increase due to a coupling of mechanical and hydrodynamic forces* 
This could be accomplished by increasing the diameter of the 
cylinder without increasing the mass very mucin This also 
would allow a coincidence of frequency at the 1st* harmonic* 

3* Take more correlation measurement in close proximity to the 
point of coincidence of frequencies to determine if there is 
any deviation from- the vs* Reynolds number trend observed 
in other regions of frequency „ 

i\ * Take correlation measurement data in the same Reynolds number 
range that o.l Earoudi did to see what changes in trends or 
values are observed for a non-rigid cylinder* 
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